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Abstract 

We extend our recent phenomenological study of pN — > NNtttt reactions to the pN — > NNtttt 
reactions for anti-proton beam momenta up to 3.0 GeV within an effective Lagrangian approach. 
The contribution of TV* (1440) with its Na decay mode is found to be dominant at the energies close 
to threshold for NNtt^tt^ and NNtt°tt° channels. At higher energies or for prnr~TT~ and NNtt^tt 
channels where A* (1440) — > Na mode cannot contribute, large contributions from double-A, 
A(1600) A*(1440)vr, A(1600) -> Avr and A(1620) Avr are found. In the near-threshold 
region, sizeable contributions from A — > Att, A — > Ntt, N — > Att and nucleon pole are also 
indicated. Although these results are similar to those for pN — > NNtttt reactions, the antinucleon- 
nucleon collisions are shown to be complementary to the nucleon-nucleon collisions and may even 
have advantages in some aspects. The PANDA/FAIR experiment is suggested to be an excellent 
place for studying the properties of relevant A^* and A* resonances. 
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I. INTRODUCTION 



As an interesting field to study baryon spectrum and properties of strong interaction, 
double pion production in pion-, photo- and electro-induced reactions has been extensively 
explored jl| . Recently, we have performed a comprehensive theoretical analysis of the double 
pion production in nucleon-nucleon collisions [2| based on the new data from CELSIUS and 

n 

COSY experiments in the past few years [3J. It is meaningful to extend the study to the 
closely related pN — > NNtttt reactions, and herein we present the results. 

The experimental studies on the pN — > NNtttt reactions were mainly performed in the 
years around 1970 Q-[lo]with some additions on the pp — > pp7i + 7i channel by the JETSET 



Collaboration at 1997 
exchange (OPE) model 



e in this channel, which can be postulated from our analysis of NN — > NNtttt 
[l2l | where the iV*(1440) was found to be important in the near-threshold region. 



JJ. The data were still scarce. On theoretical side, the one pion 
4] and Regge pole model [5], focusing on the beam momenta above 
3.0 GeV, included the double-A excitations only. However, on experimental side, there was 
an argument about the data of pp — > ppix+Tx~ at the beam momenta around 3.0 GeV [f| 
whether there was contribution from a iV* with mass about 1400 MeV and width about 
80 MeV, respectively. Also, the experiment of pp — > Pptt + tt~ at the beam momentum 
of 2.5 GeV [7J claimed an enhancement at a Air invariant mass of 1370 MeV. From the 
modern point of view, these data might show the presence of the Roper resonance iV*(1440) 
in the pp — > Pptt + tt~ channel. As a matter of fact, the iV*(1440) resonance should play 
essential ro 
reactions 

The A(1600) and A(1620) resonances are also expected to show up in the pN — > NNtttt 
reactions at high energies because they are found to be important to describe the data 
of NN — > NNtttt reactions for the beam momenta around 3.0 GeV j^J. Up to now the 
properties of these resonances are not well established and especially the nature of iV*(1440) 
is still in controversial |3_|. Therefore it is meaningful to examine whether NN — >■ NNtttt 
reactions could supply us with useful information. Also the pN — > NNtttt channels could 
serve as a complementary place to test and verify the results of pN — > NNtttt reactions. 
As we shall demonstrate later, some channels in antinucleon-nucleon collisions may be very 
suitable to settle down the problems found in nucleon-nucleon collisions. 

Our paper is organized as the following. In Sect. [TTJ we present the formalism and 
ingredients in our calculation. Then we give our numerical results and discussion in Sect. IHH 
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and a brief summary in Sect. HVl 

II. FORMALISM AND INGREDIENTS 
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FIG. 1: Feynman diagrams for NN — > NNtttt. The solid, dashed and dotted lines stand for 
the (anti)nucleon, mesons and intermediate <j(or p)-meson. The shading histograms represent 
the intermediate baryon resonances or off-shell (anti)nucleon. In the text, we use R — > NM, 
Rl -> R2M and double-^ to label (1)(2), (3)(4)(5) and (6)(7)(8), respectively. 



The Feynman diagrams of NN — > NNtttt we considered are depicted in Fig. [TJ In the 
case of iViV — > NNtttt reactions it is needed to symmetrize the initial and final nucleons so 
there are additional exchanged diagrams, which do not appear in NN — > NNtttt channels. 
The pre-emission diagrams are found to be small in NN — > NNtttt reactions Here we 
include them only for completeness. The formalism and parameters are nearly the same 
as those used in the study of the NN — > NNtttt. The commonly used Lagrangians for 
Meson- (anti)Nucleon-(anti)Nucleon couplings [13| are, 

C wNN = -^N l5lfl r-d^N, (I) 

£.aa = ^A^ 757m /C • 3"7rA„, (2) 
3 



£> v nn = -ig V NNN-f 5 r]N, (3) 
£aNN = gaNwNaN, (4) 
C pNN = -g P N N N{ln + 7T~ — (J ^vd u )T ■ p»N. (5) 

In the above and following, we explicitly specify the isospin structure of the isospin 3/2 
fields. The isospin transition operators X and /C are defined as, 



l=0,±l 



JC mn = E ( 1 h n \^ m )ei ( 7 ) 



2 '2 
z=o,±i 



where m and n are the third components of the isospin projections, and f the Pauli matrices. 
At each vertex a relevant off-shell form factor should be used and we take them as |l3| . 



/A 2 - m 2 \ n 



f m (Hf) = ( -ri — ipr J , (8) 

with n=l for 7T- and //-meson and n=2 for p- meson, fcj^, 771m and are the 4- momentum, 
mass and cut-off parameters for the exchanged meson, respectively. The coupling con- 
stants and the cutoff parameters are taken as Q: f% NN /4ir = 0.078, g^ NN /4-n- = 0.4, 
gl NN /4:TT = 5.69, g 2 pNN /iir = 0.9, = A,, = 1.0 GeV, A. = 1.3 GeV, A p = 1.6 GeV, and k 
= 6.1. We use f w AA = ^/wNn/S from the quark model The mass and width of a-meson 
are adopted as 550 MeV and 500 MeV, respectively. 

In the NN -> NNtttt reactions we have shown that the N*(1440), A(1232), A*(1600) 
and A*(1620) resonances play the major role in the considered energies [21 . Other resonances 
give negligible contributions so we can safely ignore them. The effective Lagrangians for the 
relevant resonance couplings are 

CttNA = g^NANl ■ d^nAf, + h.c, (9) 

£*NN* 1440) = g^NN^N-f^^T ■ 8^N* {im + h.C, (10) 

£ <^ ( *44o) = #^1440) NaN (im + h.c. , (11) 

ArAiV* 1440) = QvANr^AvX ' d»TtNf im + h.C, (12) 



\i440) ^ u /,JV (1440) 



ArATA* Q) = g,NA U NX ■ 7fA^ 1620) + h.C, (13) 
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£ p^* 1620) = 9 P NA* {1620 N l5 ( lfl - -^-)f • pM^ 1620) + /i.e., (14) 



^AA* 162Q) = ^AA* 162Q) \ l5 jC ■ ^7?Af 1M0) + /l.C, (15) 
£ ^a* 1600) = g*NA Uoo) Nl ■ 3"7?Af 1600)/1 + /l.C, (16) 
£ -aa* 1600) = <7*aa ?180O) A%£ • vfA^ 1600)At + /i.e., (17) 

^(Wo)A? 160 o) = toN-A^N*! ■ 0"7?A? 16OO) „ + /l.C, (18) 

For the Resonance-Nucleon-Meson vertices, form factors with the following form are used: 

^(^(i^y, (Mi 

with n=l for N* resonances and n=2 for A* resonances. We employ A* = A* = A* = A* 
= 1.0 for N*(1440), A(1232), A*(1620) and A; = 0.8 for A*(1600). The Blatt-Weisskopf 
barrier factors B{Qn*a-k) are used in the A^*(1440)-A-7r vertices [TJ 



Qn*Att + 

Here Qo is the hadron scale parameter, Q = 0.197327/i? GeV/c with R the radius of the 
centrifugal barrier in the unit of fm and chosen to be 1.5 fm to fit the data of AW — > NN7T7T 
reactions. Qn*Att and Pn*a-k is defined as, 

Qn*Ait - 1 Z S A, (21) 

JV 

2 _ K, + m|-m^ , 

" 4>< mA ' (22j 

with s x being the invariant energy squared of x particle. We introduce the Blatt-Weisskopf 

barrier factors only for A^*(1440)-A-7r vertices because other resonances, namely A* (1600) 

and A* (1620), begin to contribute at high energies so these factors have little influence on 

their behavior at the considered energies. On the other hand, as we have addressed, the 

data of nucleon-nucleon collisions at high energies are scarce so it is meaningful to decrease 

the adjustable parameters by using fewer form factors. If we would have enough data or 

go to higher energies it should be certainly necessary to include these form factors in the 

model. 

Because the mass of cr-meson is near the two-7r threshold, the following Lagrangians and 
form factor are employed for the a-ir-ir vertex 

ZaTvn = ga^d^jf ■ d^ira, (23) 



Awttt = gpirnTT X d^TT • ^, (24) 



v A 2 + q 2 j 

where q is the relative momentum of the emitted pion in the center of mass system. We use 
A = 0.8 GeV and A 2 , = 0.12 GeV 2 to normalize this form factor to unity when vr- and a- 
meson are all on-shell. The decay width of o — > tttt and p — > tttt yield g 2 n7T = 6.06 and g 2 
= 2.91. 

TABLE I: Relevant parameters used in our calculation. The masses, widths and branching ratios 



(BR) are taken from central values of PDG 
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except the BR for A* (1440) -> Avr. 



Resonance 


Pole Position BW Width Decay Mode Decay Ratio g 2 /4tt 


A*(1232)P 33 


(1210, 100) 


118 


Avr 


1.0 


19.54 


A*(1440)Pn 


(1365, 190) 


300 


Avr 


0.65 


0.51 








Act 


0.075 


3.20 








Avr 


0.135 


4.30 


A*(1600)P 33 


(1600, 300) 


350 


Avr 


0.175 


1.09 








Avr 


0.55 


59.9 








A*(1440)vr 


0.225 


289.1 


A*(1620)S 3 i 


(1600, 118) 


145 


Avr 


0.25 


0.06 








Np 


0.14 


0.37 








Avr 


0.45 


83.7 



The form factor for the baryon resonance R, Fn(q 2 ), is taken as, 

F R {q 2 ) 



A 4 



(26) 



A 4 R + (q 2 -M 2 R r 

with Ar = 1.0 GeV. The same type of form factors are also applied to the nucleon pole with 
re propagators of the exchanged meson, nucleon pole and resonances can 



A N = 0.8 GeV. 
be written as 
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G a {k a ) 



k 2 - ml + im a Y a ' 



(27) 
(28) 



Gf(M = - ^-y , (2 9 ) 

G„( 9 ) = =*^±qg) . (30) 
<T - m z N 

1/2 _ -z(^±Mg) 
G * ^ " ^MfTiM^" (31) 

-i{4±M R )G^ v {q) 
° R {q) ~ <?-Ml + iM R Y R (32) 

Here "±" is for particles and antiparticles, respectively. T R is the total width of the corre- 
sponding resonance, and G^ v (q) is defined as, 

11 2 

G^iq) = -g^u + -l^lu ± 3j^(Tm?«' ~ + Wv, (33) 

Because constant width is used in the Breit-Wigner (BW) formula, we adopt the pole posi- 
tions of various resonances for parameters appearing in the propagators. 

The coupling constants appearing in relevant resonances were determined by the empirical 
partial decay width of the resonances taken from Particle Data Group (PDG) 20j, and the 
detailed calculations of q p nr and g a NR from the R — >■ Np(a) — > Ntctt decay were given in 
Ref. [2l|. The values of cut-off in form factors were adjusted to fit the data of NN — > NN7T7T 
reactions by hand j^J. Here we would like to mention that in our fit we first determined the 
cut-off values of A and A* resonances by the pp — > nn7i + 7T + channel which had very small 
iV* contribution to be much cleaner, and then it was much easier for us to pin down the 
iV* contributions by a large amount of data in pp — > pptt + tt~ and pp — > pp7i°7r° channels. 
We tried to use the same values in the same kind of form factors for all resonances with the 
aim to reduce the number of free parameters in the model. Take the resonance form factor 
in Eq. [2(2 for example, we employed A# = 1.0 GeV for all the resonances. But in some of 
form factors, we used different cut-off values for resonances in order to reproduce the data 
better. For instance, the A* for A* (1600) in EqJT9l was not the same with other resonances. 
It should be noted that we adopted a nearly half of the decay width of A^*(1440) — > Atv 
in PDG as the recent data of NN — > N Ntctt and 7p — > ptt tt reactions favored js), H ]. 
The used decay width of A^*(1440) — > Na is the same with the value in PDG because we 
achieved an agreement with the data by adjusting the relevant cut-off parameters. So in our 
model a larger decay width of A^*(1440) — > Na compared to PDG was not required. The 
values of coupling constants and cut-off used in our computation are compiled in Table [U 
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together with the properties of the resonances and the central values of branch ratios. As we 
addressed^the parameters in Tableware the same as we used in the analysis of iViV — » NNtttt 

h 

reactions [2J. So we do not introduce any further free parameters and the calculated results 
of NN — > NNtttc in fact can be viewed as the predictions of our model. 

The amplitudes can be obtained straightforwardly by applying the Feynman rules to the 
diagrams in Fig. [TJ Isospin coefficients are considered in different isospin channels. We do 
not include the interference terms among different diagrams because their relative phases 
are unknown. The Valencia model seems to show that such terms are very small 12j, and 
our analysis of NN — > NNtttc reactions also reproduce the data well without including these 
terms. The multi-particle phase space integration weighted by the amplitude squared can 
be performed by a Monte Carlo program using the code FOWL from the CERN program 
library 23]. 



III. NUMERICAL RESULTS AND DISCUSSION 



Fig. [2] and Fig. [3] demonstrate our calculated total cross sections of four isospin channels 
together with the existing data js-ll|. Our numerical results give an overall good reproduc- 
tion to all three initial pn channels, but overestimate the initial pp channel. This may be 
due to the fact that we have not treated the initial state interactions (ISI) properly. The 
ISI usually has a weak energy dependence for the meson production processes, so adjusting 
cutoff parameters in the form factors may partly account for it effectively as for the pp 
collision j^, [if]]. However, while the pn is a pure isospin- vector state, the pp is a mixture of 
isospin-scalar and isospin-vector. The annihilation rate for the isospin-scalar NN is empir- 



ically found to be bigger than the isospin-vector by a factor a 



annihilation rates will cause different ISI reduction factors 25, 



bout 1.7 [2j. The different 



261 ] . This effect is not taken 



into account in our model calculation. The final state interactions (FSI) may also cause 
smooth energy-dependent modifications to the total cross sections [27] . Although the ISI 

-H, they 



25 



and FSI could be taken into account by some more complicated approaches 
are still of some model dependence. Since in this work we mainly investigate the relative 
importance of various resonance contributions, we have not included complicated treatments 
of ISI and FSI which are not expected to influence our main conclusions. 



In the following, we shall first address the pn — > pun n channel because it is similar to 



S 



the pp — > nmx + Tx + channel and has negligible N* contribution to be more clean. Then we 
shall discuss other three channels. We use the same definitions of various differential cross 
sections as those used in NN — > N Ntxtx reactions 

aa 

. The My and are the invariant 
mass spectra, and the angular distributions are all defined in the overall center of mass 
system. The values of vertical axis in the presented figures are all arbitrarily normalized. 
For concreteness we list the definitions of the angular distributions in the following, 

9m: the scattering angle of M; 

5if the opening angle between i and j particles; 

G-- 7 : the scattering angle of i in the rest frame of i and j with respect to the beam axis; 
the scattering angle of i in the rest frame of i and j with respect to the sum of 
momenta of i and j, corresponding to Qif defined in Ref. (jj. 

We try to give adequate information to the future measurements so we show a lot of 
observables predicted by our model in the following. PANDA is expected to install a Air 
solid angle detector with good particle identification for charged particles and photons to get 
the data of differential cross sections with good quality. Then if there are any experimental 
results in the future we can immediately know whether our model works and which aspect 
should be improved in view of the shortcomings of our model. Taking M N $ and M^ n 
as examples, we could identify the role of final state interaction and various resonances, 
respectively. The angular distributions are also useful to identify different contributions, 
especially some of which may be sensitive to the details of reaction mechanism. 



A. The channel of pn — > pmx ir 

For this channel, the N* — > Na and N* — > Atx cannot contribute. The A — > Ntx — > Ntxtx 
term is dominant for the energies below 2000 MeV in this channel as shown in Fig. &[c) and 
Fig. EJ^c). Because the A — > Ntt — > Ntxtx contribution is found to be important to describe 
the data of various NN — > N Ntxtx reactions simultaneously [2], it would be very useful to 
find other place to get some constrain on this term. The pn — > pfnx~Tx~ reaction is just 
an excellent place for such purpose. In Fig. H] we give the differential cross sections at the 
beam momentum of 1800 MeV. The peak of invariant mass spectrum M nn - is obviously 
different from that in pp — > nmx + TX + channel which makes it easy for us to identify the 
A — > Ntx Ntxtx term. The steep rise of angular distribution Q n in forward angle is 
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distinct from the symmetric shape in pp — > nmi + 7i + . This is trivial because amplitudes are 
not symmetric in the exchange of the (anti)nucleons in antinucleon-nucleon collisions. 

For the energies above 2300 MeV the A* (1600) and A* (1620) terms become significant 
and it is a good place to study the properties of them. The contribution from A* (1600) — > 
7V*(1440)7r term begins to take over as the biggest one in this energy region. 

As pointed out in our analysis of NN — > NNnn, the pp — > nnn + n + is very crucial in 
determining the cut-off parameters for the form factors of relevant A* resonances due to 
the fact that this channel has negligible N* contribution. Unfortunately, the current data 
of differential cross sections of pp — > nnn + 7r + suffer large uncertainties |3| because among 
its final four particles there are two neutrons which are difficult to detect. Especially, it is 
hard to figure out that whether there is aay damp hump structure or not to AW from 
the current data oi pp — > nn7r + 7r + [2J. The channel of pn — > pnn tt is just the same as 
clean as the pp — > nniT + i{ + but with only one antineutron in its final four particles, which 
can be easily reconstructed by the missing mass spectrum. Another ambiguity may rise up 
from the spectator proton when deuteron target is used to analyze pd — > pnp spec TT~7c~ , but 
spectator model is repeatedly confirmed to be reliable in (anti)nucleon-nucleon collisions. 
So pn — > pf?,7r~7r~ reaction is strongly suggested to be analyzed in PANDA-FAIR and it will 
be very helpful to distinguish different models. 



B. The channel of pp — > ppir + ir 

This channel is interesting because its double-A contribution mainly comes from the 
simultaneous A and A ++ excitation. As depicted in Fig. |2^a) and Fig. [3]^a), for the 
energies below 1800 MeV the iV*(1440) — > Na term gives the largest contribution while 
the nucleon pole and N — > An terms also influence the near-threshold region significantly. 
For the energies above 1800 MeV the double-A term takes over to be the most important 
one while iV*(1440) — > Na and iV*(1440) — > An rank the second and third, respectively. 
So unlike the pp — >• ppir + n~ and pp — > ppn 7r channels, in the whole energy region pp — > 
pp7T + ii~ is not suitable to extract the decay widths of iV*(1440) because of the large double- 
A contribution. However, as the best measured channel in antinucleon-nucleon collisions, 
it is useful to test models. As shown in Fig. [2]^a) and Fig. Ela), our results overestimate 
the data of the pp — > ppn + ir~ channel for the beam momenta below 2.4 GeV. As discussed 
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at the beginning of this section, this may be caused by the ISI and FSI which we have not 
included in our model calculation. The effects from ISI and FSI will not influence much our 
estimation of relative contributions from various intermediate baryon resonances. This can 
be checked by the full phase space measurement at PANDA/FAIR. 

In Fig. |S]and Fig. E]we show the calculated differential cross sections at the beam momenta 
of 1800 MeV and 2200 MeV, respectively. The iV*(1440) ->■ Na, iV*(1440) -> Att and 
double-A contributions are comparable and have important contributions. Our results are 
compatible with the old bubble chamber data measured at these energies Very similar to 
the NN collisions, the tot system is sensitive to the change of the contributions as can be seen 
in the M n +^- and cosv*™ spectrums. The double hump structure in M n + n - caused by the 
iV*(1440) — y An is obvious. The data of iViV collisions did not support these structures ^ 
but the old bubble chamber data of pp — y ppTT + 7T~ gave obvious double hump in M n + n - 
spectrums, especially at the beam momentum of about 1800 MeV 8J. Unfortunately, the 
statistics was very low and the number of selected events for each beam momentum was at 
most several hundreds, so the measured results were inconclusive. On theoretical side, the 
interference terms between iV*(1440) — > Na and iV*(1440) — > An might be relevant because 
their role on the nn invariant mass distributions have been found in nN — y nnN [28] and 
NN — y dnn [29|, so these terms should be treated with care in the future work. It should be 
mentioned that this problem may be related to the ABC effect of double-pion production in 
nuclear fusion reactions [29] , so it is meaningful to extensively study it both experimentally 
and theoretically. The luminosity of PANDA/FAIR is high enough to get the required 
production rates so the unsettled problem of the nn system can be further explored in 
pp — y ppn + n~ channel. 



C. The channel of pn — > pnn + n 

In Fig. EJ^b) and Fig. [3]^b), the iV*(1440) — y Na term is found to dominate in the whole 
considered energies and the nucleon pole term also gives significant contribution in the near- 
threshold region. It is worth to point out that unlike the pn — y pn7i + ii~ channel, the isovector 
excitation of TV* (1440) in pn —y pnn^-n^ is not enhanced compared to the pp —y pp7T + 7T~ 
because charged meson exchange is not allowed in both channels. So in a wide energy region 
it is suitable to explore the isoscalar excitation of iV*(1440). 
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The A* (1440) — > An term is the second largest for the beam momenta above 1500 
MeV and other contributions are much smaller. As can be seen in Fig. [7J the angular 
distributions of $™ and Op at 1800MeV are sensitive to the presence of TV* (1440) — > A7r 
term. Though there is possible ambiguity from the spectator proton when deuteron target 
is used in the experiment, this channel can be a better place to determine the partial decay 
widths of A* (1440) than pp — > pp7r + 7r~ and pp — > pp7r°n° where they are complicated by 
other contributions, such as the double-A and nucleon pole terms 



D. The channels of pn —> ppir tt 

The iV*(1440) — > Na term is not present in this reaction because the u-meson cannot 
decay to 7r~7r°. The double-A term is the most important one in a wide energy range as 
shown in Fig. [2^d) and Fig. [3](d). The A — > Air and A — > Ntt — > Ntttt terms have significant 
contribution below 1600 MeV and also have some contribution at higher energies together 
with the A*(1600) and A*(1620) terms. The agreement with the data is very good. 



E. The channels of pp — > ppir tt , pn — > pmr tt , pp — > npir ir 

There are no data on these three channels yet. They can be measured by PANDA ex- 
periment. The amplitudes for the pp — > ppTT°ir° and pn —>• pmr ?: channels are the same 
except for the difference of pp and pn FSI. A simultaneous measurement of these two chan- 
nels may help us to understand the pp and pn FSI. These two reactions are also similar to 
the pp — > pp7r°7i° reaction except for a different pp FSI and interference between amplitudes 
by various meson exchanges. Similarly, the pp — > npir^TT reaction has an analogous am- 
plitude with pn pp7T~TT° reaction except for the np and pp FSI and interference between 
amplitudes. A simultaneous study of these pN and pp reactions may help to pin down 
contributions from various meson exchanges. 



IV. SUMMARY 

In this paper, we present an analysis of four isospin channels of the pN — > NNnir reactions 
for the beam momenta of up to 3.0 GeV within an effective Lagrangian approach. The model 
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parameters determined from the NN — > NNnn reactions are used without introducing 
any further free parameters. We include contributions from the A*(1440), A, A*(1600), 
A* (1620) and nucleon pole to give a reasonably explanation of the measured total cross 
sections. The role of the A* (1440), A* (1600) and A* (1620) resonances in NN NNirn 
reactions have never been explored in previous studies. We give some typical differential 
cross sections which can be tested in the future measurements. We stress that PANDA (anti- 
Proton ANnihilation at DArmstadt) Collaboration at the GSI-FAIR (Facility of Antiproton 
and Ion Research) could play an important role in the baryon spectrum and a large amount 
of events on final states with baryon and antibaryon should be analyzed to extract the 
properties of relevant resonances. The conclusions reached from our model would be helpful 
to the future experiments performed at PANDA/FAIR. 
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FIG. 2: Total cross sections of NN — > NNtttt. The solid, dash-dot-dotted, dashed, dotted, 
dash-dotted, and bold solid curves correspond to contribution from double- A, A r *(1440) — > Na, 
iV*(i440) — > Att, A — > Air, A — > Ntt, and the full contributions, respectively. The data are from 



Refs. 
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FIG. 3: Total cross sections of NN — > NNtttt. The dashed, dash-dotted, dotted, dashed-dot- 
dotted, solid, and bold solid curves correspond to contribution from A* (1600) — > An, A* (1600) — > 
A^*(1440)7r, A* (1620) — > Air, nucleon pole, N — > An, and the full contributions, respectively. The 



data are from Refs. 
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FIG. 4: Differential cross sections of pn — > pnTT~ir~ at beam momentum 1800 MeV. The dashed, 
dotted and solid curves correspond to the phase space, double-A and full model distributions, 
respectively. 
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FIG. 5: Differential cross sections of pp — >• ppTr + ir~ at beam momentum 1800 MeV. The dashed, 
dotted and solid curves correspond to the phase space, TV* (1440) — s> Na and full model distribu- 
tions, respectively. 
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FIG. 6: Differential cross sections of pp — > ppTr + ir~ at beam momentum 2200 MeV. The dashed, 
dotted and solid curves correspond to the phase space, double- A and full model distributions, 
respectively. 
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FIG. 7: Differential cross sections of pn — >■ pn7r + 7r~ at beam momentum 1800 MeV. The dashed, 
dotted and solid curves correspond to the phase space, TV* (1440) — > Air and full model distribu- 
tions, respectively. 
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